Data on Neutron Widths do not refute Random— Matrix Theory 
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Recent analyses of the distribution of reduced neutron widths in the Nuclear Data Ensemble 
(NDE) and in the Pt isotopes claim strong disagreement with predictions of random-matrix theory. 
Using numerical simulations we identify the causes for the apparent disagreement. For the NDE it 
is a strong bias inherent in the method of analysis. In the Pt isotopes, it is likely to be a threshold 
state. 
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Recent analyses of the distribution of reduced neutron 
widths in the Pt isotopes and in the Nuclear Data En- 
semble (NDE) have cast serious doubt on the validity of 
random-matrix theory (RMT) in nuclei. For the Pt iso- 
topes, the analysis rejects agreement with RMT with a 
statistical significance of at least 99.997 per cent proba- 
bility For the NDE, the corresponding figure is 99.17 
per cent [1] . These results have found wide attention 0- 
|6| eroding, as they seemingly do, one of the cornerstones 
of the statistical theory of nuclear reactions 0, HJ- In 
this Letter we show that the method of analysis used in 
Ref. introduces a strong bias. The bias is responsible 
for the claimed disagreement of the NDE with RMT. In 
the Pt isotopes, a single-particle threshold state is prob- 
ably responsible for the disagreement. 

We briefly describe the method of analysis used in 
Ref. [2L 1245 reduced neutron widths in the NDE of 
Refs. H 03 

for 24 nuclides each containing between 14 
and 178 resonances were used. A cutoff procedure 11 | 
was used to eliminate p-wave resonances. The cutoff 
is based on the fact that s-wave (p-wave) resonance 
widths have an intrinsic energy dependence E 1 ! 2 (E 3 / 2 , 
respectively). Here E is the energy counted from neu- 
tron threshold, and the factors E x l 2 and E 3 ! 2 have to 
be taken at the resonance energy. The transition to re- 
duced widths removes the E 1 / 2 dependence of s-wave 
resonance widths. The remaining linear energy depen- 
dence of p-wave resonance widths was used to suppress 
these by introducing a cutoff that depends linearly on 
resonance energy and removing all widths smaller than 
the cutoff. The distribution of the remaining widths was 
analyzed using a maximum likelihood (ML) method. The 
method tests for agreement with the Porter-Thomas dis- 
tribution (PTD) predicted by RMT. The PTD is a x 2 - 
distribution with a single degree of freedom (y = 1). For 
comparison, distribution functions with other values of v 
were admitted. The test determined which of these gave 
best agreement with the data. This led to the above- 
mentioned rejection of RMT. 



The normalized ^-distribution for v degrees of free- 
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Here (x) denotes the average of x, and V is the complete 
gamma function. When the distribution is cut off at 
Xmin > the resulting normalized function g cut has the 
form .g cu t(x, v, (x)) = (1/C)g(x, v, (x))0(x - x min ) where 
6 is the Heaviside function. With T — x m i n /(x) the 
normalization constant C is given by 



C = 
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where T(x,y) is the incomplete gamma function. Sup- 
pose now that for some fixed values of v and (x), 
9cut{x,v, {x)) is the correct distribution function for the 
variable x after cutoff. Then the normalized probabil- 
ity for x to have the measured value Xj at a resonance 
with energy Ei is (1/C,)#(xj, v, (x))8(x,-x min (i)). Here 
x m in(i) is the value of the cutoff at Ei, and Cj is given by 
Eq. ([2]), with T — x m i n (i)/(x). For a set of iVi nc i reduced 
widths Xi of resonances with energies Ei, i = 1, . . . , Ai nc i 
remaining after the cutoff (so that 9(xi — x m i n (i)) = 1 for 
all i = 1, . . . , Aind) and included in the ML analysis, the 
likelihood function L is the product of the probabilities 
for finding resonance width Xj at energy Ei, 



(3) 
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Now the argument is turned around: One considers L 
as a function of the two variables v and (x) and argues 
that the maximum of L determines the most likely val- 
ues of these parameters. The functional dependence of 
L on v and (x) in the vicinity of the maximum deter- 
mines the error and, thus, the statistical significance with 
which a value of v that differs from the value at the max- 
imum, is rejected. That is how the figure of 99.17 per 
cent cited above for the NDE is arrived at. In actual 
fact the method used in Ref. @ is slightly different. In- 
hdstej^lietitteh^l^hMQMiaMiation factors Ci, a single joint 
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normalization factor was used. Wc return to that point 
below. 

We have tested the method by simulation. Using a 
random-number generator we have drawn N real num- 
bers yi, i — 1,...,N from a Gaussian probability dis- 
tribution with unit width. The squares determine N 
widths Xi = yf. In the limit N — > oo, these by con- 
struction obey the PTD with (x) = 1 . (We have also 
simulated widths for (x) — 0.1 and for (x) = 10, with 
no discernible differences in the analysis.) Each of these 
widths is associated with one of N resonance energies 
Ei = id with constant level spacing d [l2j]. We repeat 
that procedure 2500 times, thereby generating for each 
value of N an ensemble of widths and associated ener- 
gies. For every member of the ensemble, we use the cutoff 
and apply the ML analysis described above. For clarity 
we distinguish the input value (x) = 1 used to generate 
the data from the result of the analysis. The maxima 
of the 2500 likelihood functions L are located at values 
of v denoted by Vk and of the average width denoted by 
(r°)fc, k = 1, . . . , 2500. The uncertainty of was es- 
timated by determining those values of v where In L has 
decreased by the factor 1/2. Since L is slightly asymmet- 
ric, <jfe was taken as the average of the two value s. We 
calculated the mean values V of the and (r°) of the 
widths. We determined the statistical error er cns of the 
ensemble average V from a~^ s = y\ ctj7 2 . 

A cutoff on the widths is defined by elimination of all 
normalized widths that are smaller than the cutoff value 
^min- The normalized widths are the actual simulated 
widths Xi divided by their average. For finite N that 
average differs slightly from unity. The ML analysis is 
then applied to the remaining (unnormalized) iVnci < N 
widths and the associated remaining energies. The latter 
do not have constant spacing. The likelihood function is 
given by Eq. (|3]), with (x) replaced by the unknown width 
(T°), with C t given by Eq. ©, and with T = x miu /(r° n ). 
We have used six different cutoff procedures, (i) No cut- 
off (a; m i n = 0). (ii) A constant cutoff cc m i n independent 
of resonance energy as in Ref. [11]. Then all Cj have 
the same value ^j. (iii) A linear cutoff x m [ n (Ei) = 
X m i n Ei/E max where E max = Nd and where X min mea- 
sures the slope. With T(Ei) — x min (i?j)/(r°) the nor- 
malization constants are Ci = T(y/2, vT{Ei)/2)/T(v/2). 
(iv) A linear cutoff as under (iii) but for a modified form 
of the likelihood function ^ where every C, is replaced 
by C, the linear smooth average of the Ci, 
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of L actually 

used in Refs. (EH. (v) A cutoff given by x m i n (Ei) = 
X min [Ei/E max \ 3 / 2 . That is the cutoff used in Ref. 
for the reduced neutron widths in the Pt isotopes. The 
normalization constants Ci have the same form as under 
(iii). (vi) Same as under (v) but with every Ci replaced 
by C as defined in Eq. (fj|. Method (vi) was used in 
Ref. [l[. Method (iii) (method (v)) corresponds to the 



standard ML approach and is supported by mathematical 
theory while method (iv) (method (vi), respectively) is 
not. We test whether methods (iv) or (vi) yield results 
that differ from those of methods (iii) or (v) , respectively. 

We have simulated data for N = 100, N = 450, N = 
1100. Here N = 100 corresponds roughly to a typical set 
of widths for a given nuclide in the NDE, N = 1100 is 
approximately the total number of widths in the NDE, 
and N — 450 provides an intermediate value. The three 
values of TV are also used to test whether there is a sys- 
tematic TV-dependent trend in the ML analysis. Some 
of our results are shown in the tables. In Table I we 
show for methods (i) to (iv) and the two extreme values 
of N the cutoff parameters as defined above, the average 
number 7V; nc i of widths retained after cutoff, and the en- 
semble averages V and (F°) with their statistical errors 
resulting from ensemble averaging. Since we used 2500 
simulations, the errors are small in all cases. Method 
(i) serves as a test of our method. Method (i) yields for 
all three 7V-values V m 1 and (r°) w 1, consistent with 
the input. The values of V and (r°) determined with 
methods (ii) to (iv) increase with increasing cutoff and 
decrease very slowly with increasing N. Methods (iii) 
and (iv) give very similar results for all three values of 
N; the symmetrization Q does not adversely affect the 
results. In view of the small statistical error, agreement 
with RMT would be firmly excluded for all non-zero cut- 
offs that are bigger than 0.02 - in spite of the fact that 
the data have been generated from the PTD. These facts 
establish the existence of a strong cutoff-dependent bias 
in the ML analysis of neutron resonance widths [3] . 

How does the bias arise? How can it be that the PTD 
used to generate the data does not give the best fit to 
these same data? Let us consider method (ii) (constant 
cutoff at a; m in). An unbiased analysis would test whether 
after removal of all widths smaller than a; m ; n , the shape of 
the distribution of the remaining widths is that of a PTD 
with a cutoff at x m - ln . The analysis might, for instance, 
test for maximum overlap, irrespective of normalization. 
The ML analysis, in contradistinction, assigns a normal- 
ized probability to each width and searches for the max- 
imum of the product of these probabilities. That intro- 
duces the factor C^^ 1 = [T(y/2, vT /2) /T{v /2)]- N ^ 
in Eq. ([3]). That factor has a non-trivial dependence on 
v and, through T, on (r°). The maximum of L is, thus, 
not determined by maximum overlap. It is influenced by 
the normalization factors. That causes the bias. Similar 
considerations apply for non-constant cutoffs (methods 
(iii) to (vi)). 

The statistical error on V is small and would decrease 
further if the number of simulations were increased. How- 
ever, the maximum of L for an individual simulation 
with index k does not, in general, coincide with V. For 
the cases considered in Table I, Table II shows the root- 
mean-square deviations 5v of the Vk from V. Such devia- 
tions are likewise expected to occur for a set of N widths 
in any given nucleus. Table II also shows the average a 
of the widths Ofc defined above. We note that a > 5v in 
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TABLE I: Ensemble averages v and (r°) for different sample sizes, different methods, and different cutoff parameters. 









Method (ii) 




Method (iii) 


Method (iv) a 


N 


^rain/ ^"min 


iVincl C 




\ l nl 


M. , ac 

1 Mild 


v <r°> 


v (r°) 






100 


n QQ4 4- n nri9 

u.yyi zn u.uu^ 


u.yuo zn u.uuo 




Same as Method (ii) 


Same as Method (ii) 




0.01 


92 


n QQ7 4- n nru 

u.yyi zn u.uu*± 


i.u^o zn u.uuo 


95 


0.999 ± 0.004 1.016 ±0.003 


1.002 ± 0.004 1.016 ± 0.003 




0.02 


89 


i niTH 4 n no4 

i.uuo zn u.uu*± 


i 07Q 4 n oo^ 

i.ui y zn u.uuo 


93 


0.999 ± 0.004 1.039 ± 0.003 


1.003 ± 0.004 1.039 ± 0.003 


100 


0.05 


83 


i n^d 4 n nn'i 

i.uoi zn u.uuo 


J-.-LOO zn u.uui 


88 


1.003 ± 0.004 1.087 ± 0.003 


1.009 ± 0.004 1.087 ± 0.003 




0.10 


75 


i ill 4 n nnfi 


1 9 C >8 4 004 

l.idO Zn U.UU4: 


83 


1.024 ±0.005 1.146 ±0.004 


1.031 ±0.005 1.146 ±0.004 




0.20 


66 


i ^o^ 4 n fins 

r.ouo zn u.uuo 


1 41 Q 4 004 
i.'iiy zn u.uu^fc 


77 


1.081 ±0.006 1.236 ±0.004 


1.089 ± 0.006 1.236 ± 0.004 




0.30 


58 


1.532 ± 0.009 


1.560 4 0.005 


72 


1.141 ±0.006 1.311 ±0.004 


1.149 ±0.006 1.310 ±0.004 




Q d 


1100 


0.9999 ± 0.0007 


0.9965 =b 0.0008 




Same as Method (ii) 


Same as Method (ii) 




0.01 


1012 


1.009 4 0.001 


1.082 ±0.001 


1042 


1.003 ±0.001 1.052 ±0.001 


1.003 ±0.001 1.052 ±0.001 




0.02 


977 


1.018 4 0.001 


1.121 ±0.001 


1017 


1.007 ±0.001 1.077 ±0.001 


1.007 ±0.001 1.077 ±0.001 


1100 


0.05 


906 


1.056 ± 0.002 


1.206 ±0.001 


971 


1.018 ± 0.001 1.128 ±0.001 


1.018 ± 0.001 1.128 ±0.001 




0.10 


828 


1.137 4 0.002 


1.313 ±0.001 


917 


1.041 ±0.002 1.190 ±0.001 


1.041 ±0.002 1.190 ±0.001 




0.20 


721 


1.327 4 0.002 


1.485 ± 0.001 


844 


1.100 ± 0.002 1.287 ±0.001 


1.099 ±0.002 1.287 ±0.001 




0.30 


643 


1.548 ± 0.003 


1.635 ±0.002 


789 


1.162 ±0.002 1.365 ±0.001 


1.159 ±0.002 1.365 ±0.001 



"iVjaci is always the same for Methods (iii) and (iv). 
'For method (ii), £ m i n is the relevant cutoff parameter; for meth- 
ods (iii) and (iv), X m i n is the quantity that must be specified. 
Uncertainties in A^ nc ] are much less than 1. 
%mm = or X m i n = corresponds to method (i). 



TABLE II: Error estimates for N = 100 and N = 1100 for 



individual spectra 





3>min / 


Method (ii) 


Methods (iii)-(iv) 


N 


V" a 


5v 


a 


5v 


a 




o 6 


0.13 


0.12 


Same as 


Method (ii) 




0.01 


0.20 


0.20 


0.15 


0.18 




0.02 


0.21 


0.23 


0.15 


0.20 


100 


0.05 


0.24 


0.27 


0.16 


0.22 




0.10 


0.25 


0.32 


0.17 


0.25 




0.20 


0.28 


0.40 


0.19 


0.29 




0.30 


0.31 


0.47 


0.20 


0.32 




o 6 


0.03 


0.03 


Same as 


Method (ii) 




0.01 


0.06 


0.06 


0.05 


0.05 




0.02 


0.06 


0.07 


0.06 


0.06 


1100 


0.05 


0.07 


0.08 


0.06 


0.07 




0.10 


0.07 


0.09 


0.07 


0.07 




0.20 


0.08 


0.11 


0.07 


0.08 




0.30 


0.09 


0.13 


0.08 


0.09 



"For method (ii), ic m ; n is the relevant cutoff parameter; for meth- 
ods (iii) and (iv), X m i n is the quantity that must be specified. 

'^'min = or X m ; n = is the same for each method; these are 
the results when method (i) is utilized. 



all cases where a cutoff is utilized. If all Ok were equal, 
the statistical error on V in Table I would simply be given 
by a/VK with K = 2500. 

In Ref. the weighted average of i>u over 24 nuclides is 



given as Fkoo = 1.217±0.092. Because of the small num- 
ber of nuclides and because of the spread of the number of 
resonances per nuclide, the error on Fk oc is much larger 
than for the 2500 simulations in our Table I. Within that 
error, the result for Fkog agrees, however, with the re- 
sults of our Table I for cutoff values bigger than 0.2. In 
other words, the NDE behaves exactly as is expected for 
data that follow the PTD. Thus, the conclusions drawn 
in R ef. m on the NDE cannot be upheld. The authors of 
Ref. [101 bave concluded that the distribution of neutron 
widths in the NDE agrees with the PTD. There is no 
reason to doubt that conclusion. 

The Pt isotopes are special. While the analysis of the 
NDE in Ref. Q yields for all nuclides values for that 
are bigger than unity or, within errors, are consistent 
with unity, the values for v exp given in Table I of Ref. [l[ 
for the three Pt isotopes range from 0.47 to 0.60 with 
errors not larger than 0.28. While in Ref. @ the NDE is 
claimed to disagree with RMT because FK oe is too big, 
the disagreement claimed in Ref. [l| comes about because 
u exp is too small. That suggests that the distribution of 
widths in the Pt isotopes is significantly different from 
that of the nuclides in the NDE, and that a special mech- 
anism may here be at work. The investigation of the 
Pt isotopes in Ref. [2] was motivated by the fact that 
for these mass numbers, the neutron strength function 
is close to its maximum, so that p-wave contributions 
are effectively suppressed by about a factor 10 The 
enhancement of the s-wave strength function is due to 
a threshold state: With increasing mass number, the 4s 
single-particle state of the shell model is pulled into the 
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TABLE III: Values of V for spectra with N = 161 and £ max = 
15.93 keV for different values of \Eo\ when method (vi) is 
applied. 





0.5 


\E \ (keV) 
5 


20 





0.797 ±0.001 


0.890 ± 0.002 


0.970 ± 0.002 


0.01 


0.690 ± 0.002 


0.851 ±0.002 


0.976 ± 0.003 


0.02 


0.676 ± 0.002 


0.846 ± 0.003 


0.979 ± 0.003 


0.05 


0.666 ± 0.002 


0.850 ± 0.003 


0.992 ± 0.003 


0.10 


0.674 ± 0.003 


0.866 ± 0.003 


1.011 ±0.004 


0.20 


0.704 ± 0.003 


0.908 ± 0.004 


1.057 ± 0.004 


0.30 


0.733 ± 0.003 


0.946 ± 0.004 


1.098 ± 0.004 



nuclear potential well, and it passes threshold at or near 
the masses of the Pt isotopes. In Ref. [1] it was shown 
that a threshold state at energy \E \ either above or be- 
low threshold modifies the s-wave neutron widths by the 
factor 



E 



E+ En 



(5) 



Here v E is the usual s-wave penetration factor, and 
1/(E ± |-Eb|) is the resonance enhancement due to the 
threshold state. For each resonance, f 2 (E) has to be 
taken at the resonance energy Ei. In Ref. [l[, the re- 
duced widths of the Pt isotopes were analyzed using 
an energy-dependent cutoff proportional to E 3 / 2 with- 
out taking into account the denominator in Eq. ([5]). To 
test for the influence of the threshold state, we simulated 
2500 sets of N widths Xi and N equally spaced ener- 
gies Ei in the interval {0, £ max } as described above, with 
N = 153 and E max = 4.98 keV for 192 Pt and N = 161 
and E max = 15.93 keV for 194 Pt, plus some other cases 
for additional information. For each realization of the en- 
semble, we multiplied every Xi by l/(Ei ± \Eq\). (Other 
energy-independent factors contributing to f 2 were ig- 
nored since a rescaling of all widths does not affect the 
value of Vk)- The resulting set of modified widths and as- 
sociated resonance energies was subject to the ML anal- 
ysis with an energy-dependent cutoff. 

Some of our results are shown in Table III. Without 



cutoff, significant deviations from the input value v = 1 
appear only when \E \/ '-E max < 1/2 or so. For the values 
of | So | shown, these deviations are considerably smaller 
than the experimental ones in the Pt isotopes. Methods 
(iii) and (iv) again give very similar results; for some com- 
bination of values of \Eq\ and cutoff parameter, values of 
V as small as 0.7 are found. The results of methods (v) 
and (vi) differ, method (v) generally giving smaller values 
for v. In both cases these values first drop and then in- 
crease with increasing cutoff, the values for method (vi) 
attaining a minimum around 0.70 for \Eq\ = 0.5 keV. 
That is close to but larger than the experimental values 
that range from 0.5 to 0.6. However, the replacement 
Xi — > f 2 (Ei)xi is only approximately valid if \Eq\ is as 
small as 1 keV or so. It is shown in Ref. 4] that in such a 
case, the entire i?-matrix analysis of the scattering data 
must be redone, taking the existence of the threshold 
state into account from the outset. It is not known how 
exactly that affects the distribution of the x^. We con- 
clude that a single-particle state close to threshold is a 
likely candidate for the deviations found in the Pt iso- 
topes. Establishing its existence would require, however, 
a complete re-analysis of the scattering data that would 
not make use of standard i?-matrix codes. 

In summary, we have shown that the ML analysis used 
in Ref. @ produces a strong bias. The bias distorts 
the value of V determined from the data and leads to 
an apparent severe disagreement with RMT predictions. 
When the bias is taken into account, the distribution 
of neutron widths in the NDE is consistent with RMT. 
For the Pt isotopes we have shown that a single-particle 
threshold state is a likely cause for the observed disagree- 
ment with RMT. To firmly establish this assertion, a 
complete re-analysis of the scattering data would be re- 
quired, however. 
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